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ABSTRACT 

High-magnification gravitational microlensing events provide an important channel of detecting planetary 
systems with multiple giants located at their birth places. In order to investigate the potential existence of 
additional planets, we reanalyze the light curves of the eight high-magnification microlensing events for each 
of which a single planet was previously detected. The analyzed events include OGLE-2005-BLG-071, OGLE- 
2005-BLG-169, MOA-2007-BLG-400, MOA-2008-BLG-310, MOA-2009-BLG-319, MOA-2009-BLG-387, 
MOA-2010-BLG-477, and MOA-201 l-BLG-293. We find that including an additional planet improves fits 
with Ax 2 < 80 for seven out of eight analyzed events. For MOA-2009-BLG-319, the improvement is relatively 
big with A% 2 ~ 143. From inspection of the fits, we find that the improvement of the fits is attributed to 
systematics in data. Although no clear evidence of additional planets is found, it is still possible to constrain 
the existence of additional planets in the parameter space. For this purpose, we construct exclusion diagrams 
showing the confidence levels excluding the existence of an additional planet as a function of its separation and 
mass ratio. We also present the exclusion ranges of additional planets with 90% confidence level for Jupiter, 

Saturn, and Uranus-mass planets. 

Subject headings: gravitational lensing: micro - planetary systems 


1. INTRODUCTION 

Since the first disc overies (IWolszczan & FraiTI Il992h 
iMavor & QueloJ 119951) . the number of known extra-solar 
planets has been rapidly increasing and it reaches ~ 2000 
(http://exoplanet.eu: ISchneider et al.l 1201 lb . A significant 
fraction of these planets are members of multi-planet systems 
dRowe et alJl2014t iLissauer et ~aPl20l2h . which were mostly 
discovered by the radial-velocity and transit methods. 

Multi-planet systems can also be detected by using the mi¬ 
crolensing method. Gravitational microlensing occurs by the 
chance alignment along the line of sight toward a background 
star (source) and a foreground object (lens), resulting in the 
magnification of the source brightness due to the bending of 
light by the gravity of the lens. When a lensing object is 
a star accompanied by a planet, the planet causes astigma¬ 
tism in the bending of light and induces formation of caustics 
which refers to the envelope of light rays on the source plane 
at which the magnification of a point source becomes infinity. 
A planetary lensing signal occurs when the source approaches 
caustics and it appears as a short term anomaly superposed 
on t he smooth and symmetric lensing light curve of the host 
star (iMao & PaczvnskJll99ll : lGould & Loebll 19921) . If a plan¬ 
etary system is composed of multiple planets, each planet 
induces its own caustics and the multiplicity of the plan¬ 
ets can be identified if the source trajectory pmsses multiple 
caustics ind uced by the individual planets (lHan et all 12001b 
lHan 1120051) . The sensitivity to multiple planets for general 
microlensing events is low because planet-induced caustics 
are small and thus the geometrical probability for the source 
to pass multiple caustics is low. However, the sensitivity is 
high for a subset of lensing events where the brightness of a 
source star is greatly magnified (iGaudi et al.l 119981) . This is 
because one of the caustics induced by each planet is always 
located very close to the host star around which the source 
of a high-magnification event passes by dGriest & Safizadehi 


1998). Ind eed, two multi-planet systems ide ntified by mi¬ 
crolensing (IGaudi et al.l 120081 : lHan et all 120131) were discov¬ 
ered through this channel. 

Microlensing discoveries of multi-planet systems are im¬ 
portant for better understanding of planet-formation mech¬ 
anism. According to the core accretion theory ( Ida & Linl 
i2004j) . giant planets are believed to form in group in the re¬ 
gion beyond the snow line where the lower temperature in this 
region makes many more solid grains available for accretion 
into planetesimals. However, most known planets in multiple 
systems reside well within the snow line. These close-in plan¬ 
ets are believed to h ave migrated from the outer region where 
they were formed (ID’ Angelo & Lubowl 120081: [Luhow & Ida 
1201 ll) . Detecting planets at their birth places by using the 
transit and the radial-velocity methods is difficult due to their 
low sensitivity to wide-separation planets. By contrast, mi¬ 
crolensing is sensitive to planets located around and beyond 
the snow line and thus provides a channel to study planets 
where they formed. 

In this work, we reanalyze the light curves of the eight high- 
magnification microlensing events for each of which a single 
planet was previously detected. From the analyses, we inves¬ 
tigate the existence of additional planets. 


2. SAMPLE 

In published articles, there exist 29 planetary microlens¬ 
ing events from which 31 planets were found. Among them, 
events for our analyses are selected based on the follow¬ 
ing criterion. First, we choose planetary events with high 
magnifications due to their high sensitivity to multiple plan¬ 
ets. With an adopted threshold magnification A± ~ 50, 
there exist 15 events. Among the planetary systems found 
from these events, two are already known to be in multi- 










































2 


Constraints on Additional Planets in Known Planetary Systems 


TABLE 1 

List of sample events and the properties of the previously dicovered planets 


Event 

Amax 

Previously dicovered planet 

Reference 



Mass 

Separation (AU) 


OGLE-2005-BLG-071 

59 

3-8 tt Mi 

3 6+ 0 - 2 

2 7+¥:7 
-1.4 

Udalski et al. (2005), Dong et al. (2009) 

OGLE-2005-BLG-169 

880 

13?f M® 

Gould et al. (2006) 

MOA-2007-BLG-400 

629 

0.83^3J Mj 
74f[y M® 

q 79 + 0.38 / 6 5 +3 - 2 
u '' -0.16 ' 

Dons et al. (2009) 

MOA-2008-BLG-310 

366 

1 9 C+ 0.10 
A ’ ZJ - 0.10 

Janczak et al. (2010) 

MOA-2009-BLG-319 

208 

50^ M® 

2 4+ 1 -* 
z -^- 0.6 

Mivake et al. (2011) 

MOA-2009-BLG-387 

49 

2.6ft* Mi 

1 O+0.9 
^ ‘ -0.7 

Batista et al. (2011) 

MOA-2010-BLG-477 

MOA-201 l-BLG-293 

599 

418 

1 -5 + S | Mj 
2.4 Mj 

2.0+ 30 

1 0 + 0 ' 1 / 3 4+^ 

Bachelet et al. (2012) 

Yee et al. (2012) 


planet system^ and another two are in stellar binary sys¬ 
tem s|J] We exclude the events with already-known multi¬ 
planet systems and circum-stellar binary planets from our 
analysis. Second, we choose events with good coverage of 
the peak region. This is because signatures of multiple plan¬ 
ets for high-magnification events show up near the peak of 
a light curve and thus good coverage of the peak region 
is essential to constrain additional planets. Third, we con¬ 
fine analysis to events for which there is no ambiguity in 
the interpretation of the known planetary signals. Applying 
these criteria leave 8 events including OGLE-2005-BLG-071, 
OGLE-2005-BLG-169, MOA-2007-BLG-400, MOA-2008- 
BLG-310, MOA-2009-BLG-387, MOA-2010-BLG-477, and 
MOA-201 l-BLG-293. 

In Table[I] we list the lensing events that we analyzed. Also 
presented are the peak magnifications A max of the individual 
events,the physical parameters (mass and separation from the 
host) of the previously discovered planets, and the references 
of the original analyses. In Figure Q] we present the light 
curves of the events. We note that the same data sets as those 
used in the original analysis are used in our analysis for the 
consistency of results. 

3. ANALYSIS METHOD 

In order to investigate the existence of an additional planet 
for a lensing event with an already known planet, it is required 
to conduct triple-lens modeling. For a triple-lens system, the 
image positions for given locations of the lens and source are 
found by solving the equation of lens mapping (lens equation) 
that is expressed as 






, Z~Zh,k 
=1 


( 1 ) 


Here k = 1, 2, and 3 denote the individual lens components, 
£ = C + if1, Zh,k = *l ,k + iy~L,k and z = x+iy are the complex no¬ 
tations of the source, lens, and image positions, respectively, 


1 The first known multiple microlensing planetary system is OGLE-2006- 
BLG-109L, where a host star with a mass ~ 0.5 Mq is orbited by two planets 
with masses ~ 0.71 Mj and ~ 0.27 Mj and separations from the host ^ 2.3 
AU and ~ 4.6 AU, respectively iGaudi et all2 008: Benne ttet all2010h . The 
second system is OGLE-2012-BLG-0026L, where a host star with a mass 
~ 0.82 Mq is orbited by two planets with masses ~ 0. 11 Mj and ^ 0.6 7 Mj 
and separations ~ 3.8 AU and 4.6 AU, respectively JHan et aI]l2013D . 

■t OGLE-2013-BLG-0341Lb is the first microlensing planet in a binary 
stellar system, where a terrestrial planet (~ 2 Mq) lies at ~ 0.8 AU from its 
host with a mass ~0.10-0.15 Mq and the host itsel f orbits another sta r with 
a mass ~ 0.15 Mq and a separation ~10-15 AU JGould et alj|2014l) . The 
second microlensing planet in binary is OGLE-2008-BLG-092Lb, where a 
~ 4 Mq planet orbits a ~ 0. 7 Mq star at ^ 18 A U. The host has a low-mass 
companion with ~ 0.15 Mq (Poleski et all 20141) . 


z denotes the complex conjugate of z, and ej = Mk/M tot is the 
mass fraction of each lens component (IWittll 199~0l) . Once the 
image locations are found, the lensing magnification is com¬ 
puted by 


: 57 ■Ad ■ 


acacy 1 

dzdz) z=Zi ’ 


( 2 ) 


where A, is the magnification of each image. 

Due to non-linearity, the lens equation for a triple lens sys¬ 
tem cannot be analytically solved. However, lensing magni¬ 
fications can be numerically computed by using the method 
known as the inverse ray-sh ooting method (Kav ser et alJ 
11986t[Schneider & Weisslll986i) . In this method, a large num¬ 
ber of rays are uniformly shot from the image plane, bent by 
the lens equation, and land on the source plane. Then, lensing 
magnifications are computed as the ratio of the number den¬ 
sity of rays on the source plane to that on the image plane. 
The lens equation can also be solved semi-analytically be¬ 
cause the equation is expressed as a tenth-order polynomial 
in z and the image positions and subsequent lensin g magnifi¬ 
cations are comp uted by solving the polynomial ( Rhiel 120021 : 
ISong et al.l2014l) . The semi-analytic method has an advantage 
of fast computation but it is difficult to be used for computing 
lensing magnifications affected by finite-source effects. On 
the other hand, the numerical method can be used to compute 
finite-source magnifications but computation is slow. 

Describing the light curve of a triple-lensing event requires 
ten basic parameters. The first three of these parameters de¬ 
scribe how the source approaches the lens. They are the time 
of the closest source approach to a reference position of the 
lens system, to, the normalized projected separation between 
the source and the lens reference position at that moment, uq, 
and the Einstein time scale, fg. that is defined as the time re¬ 
quired for the source to cross the angular Einstein radius of 
the lens Oy. Another six parameters characterize the lens sys¬ 
tem. They are the projected separations between the host and 
the individual planet companions, s\ and ,s' 2 , and the mass ra¬ 
tios of the companions to the host, <71 and q 2 , the angle be¬ 
tween the source trajectory and the axis connecting the host 
and the first planet, a, and the orientation angle of the second 
planet with respect to the host-first planet axis, ip. We note 
that length scales describing lensing phenomenon are usually 
expressed in units of 8e and uq, .v 1 and si are normalized by 
0 e- The last parameter is the angular size of the source 0* 
in unit of the Einstein radius, p t = 0*/0 e (normalized source 
radius). This parameter is needed to describe light curve de¬ 
viations affected by finite-source effects. 

Besides the basic lensing parameters, precise description of 
lensing light curves often requires additional parameters to 
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Fig. 1.— Light curves of planetary microlensing events reanalyzed in this work. Grey points with error bars are the data used in our analysis. Also presented 
are the model light curves base on single-planet models. 


account for higher-order effects. One such effect is caused 
by the positional change of an observer induced by the orbital 
motion of the Earth around the Sun. Considering this parallax 
effect ( iGould 1992 ) requires two additional parameters ttu,n 
and tte,e that represent the two components of the lens paral¬ 
lax vector 7 r E projected on the sky in the north and east equa¬ 
torial coordinates, respectively. Similarly, lensing light curves 
can also be affected by the change of the lens position caused 
by the orbital motion (IShin et al Jl20TTl 1201 2t ISkowron et al.1 

l20Tll: IPark et~ani20Tllh 

With the lensing parameters, we conduct triple-lens mod¬ 
eling of the events that were previously analyzed based on 
binary-lens modeling. For all analyzed events, the anomaly 
is dominated by the signal of the already reported planet 
and thus the signal of the potential additional planet would 
be small. We, therefore, set the lensing parameters related 
to the known planet (si,^i) fixed as outlined in iKubas et all 
(120081) . The existence of an additional planet is then in¬ 
spected by checking whether the two-planet model improves 
the fit with respect to the single-planet model. Not know¬ 
ing the characteristics of the second planet, we search for the 
second-planet parameters .s '2 and qi by inspecting solutions in 
wide ranges of the parameter space spanning —1 < log,s '2 < 1 
and -6 < logf /2 < 2, respectively. For a given set of (so-c/ih 
the other lensing parameters are searched for by minimizing 
X 2 in the parameter space using Markov Chain Monte Carlo 
(MCMC) algorithm. 

For computing triple-lensing magnifications, we apply both 
the semi-analytic and numerical methods. For the region near 
the peak of the light curve, we use the numerical inverse ray¬ 
shooting method because this region is likely to be affected 
by finite-source effects during the source star’s approach close 
to the planet-induced caustics. For other regions, we use the 
semi-analytic method for fast computation. Computing finite- 
source magnifications is further accelerated by applying the 


“map making” method dDong et al.ll2006i) . In this method, a 
map of rays for a given set of (si,^) and (* 2 , # 2 ) is constructed 
based on the inverse ray shooting method and then it is used 
to produce many different light curves resulting from various 
source trajectories without producing extra maps. In comput¬ 
ing finite-source magnifications, we consider limb-darkening 
variation of the source star surface brightness by modeling the 
surface brightness profile as S\ oc 1 -T A (1 — 1.5 cos <fi), where 
T A is the limb-darkening coefficients corresponding to spe¬ 
cific passband A and (j> is the angle between the line of sight 
toward the source star and the normal to the source surface. 
We use the same limb-darkening coefficients that were used 
in the previous analyses where first planets were reported. 

If higher-order effects were reported in the previous anal¬ 
yses, we also consider them. Parallax effects were reported 
for events OGLE-2005-BLG-071, MOA-2009-BLG-387, and 
MOA-2010-BLG-477. On the other hand, obvious effect of 
lens-orbital motion was reported for none of the events. We 
note that the effect of an additional planet is confined to the 
narrow region of the light curve peak while both parallax and 
lens-orbital effects influence the overall shape of a lensing 
light curve. Therefore, the existence of the second planet does 
not affect the lensing parameters of higher-order effects. 

4. RESULT 

4.1. Existence of Additional Planets 

Table[2]shows the result of the analysis presented as x 2 val¬ 
ues of the binary and triple lens models. From the comparison 
of the x 2 values, it is found that including an additional planet 
improves fits with Ax 2 < 80 for seven out of the eight ana¬ 
lyzed events. For MOA-2009-BLG-319, the improvement is 
relatively big with Ax 2 = 142.7. 

Although the improvement of the fit for the event MOA- 
2009-BLG-319 is formally significant, one cannot exclude the 
possibility of systematics in data, which is often masqueraded 
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TABLE 2 

COMPARISION OF BINARY AND TRIPLE LENS FITS AND EXCLUSION RANGE OF A 
SECOND PLANET 


Event 

Xbi 

V*. 

/ Mri 

1 * 

<1 

Exclusion range (AU) 

Jupiter Saturn Uranus 

OGLE-2005-BLG-071 

1305.6 

1283.7 

21.9 

1.5 - 5.1 

_ 

- 

OGLE-2005-BLG-169 

533.4 

518.5 

14.9 

0.6- 11.3 

1.0 - 7.1 

- 

MOA-2007-BLG-400 

770.8 

760.0 

10.8 

0.8-4.5 

1 .2-2.6 

- 

MOA-2008-BLG-310 

3154.4 

3076.4 

78.0 

0.7 - 2.0 

1.1 - 1.3 

- 

MOA-2009-BLG-319 

6180.2 

6037.5 

142.7 

0.5-7.9 

0.9-5.0 

1.6-2.9 

MOA-2009-BLG-387 

1298.4 

1269.3 

29.1 

0.9-4.4 

1.6-2.5 

- 

MOA-2010-BLG-477 

4444.0 

4394.4 

49.6 

0 .8-12.1 

1.7-5.8 

- 

MOA-201 l-BLG-293 

5089.2 

5069.9 

19.3 

1.0-3.5 

- 

- 


NOTE. — The \ 2 difference represents A\ 2 = Xbi — where and xjn are for the 
binary and triple lens models, respectively. The exclusion ranges are estimated based on 
90% confidence level 
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Fig. 2.— Comparison of the triple and binary models of MOA-2009-BLG- 
319. The lower two panels show the residuals from the individual models. 



as planetary signals. In order to check the genuineness of 
the additional planetary signal, we inspect the fits from the 
triple and binary lens models. Figure [2] shows the peak por¬ 
tion of the light curve. We plot data sets from MOA, CTIO, 
and FTN observations for which most signal of the additional 
planet comes from: A% 2 = 78.4 for MOA, A% 2 = 11.0 for 
CTIO, and A% 2 = 35.5 for FTN data. Also plotted are the light 
curves of the triple and binary models. It is found that the two 
models cannot be distinguished within the resolution of the 
plot. To better show the difference between the two models, 
we also present the residuals from the individual models. We 
find that the improvement by the triple-lens fit is ~ 0.01 mag 
level which is equivalent to the scatter of the data, implying 
that the signal can centainly be attributed to the systematics in 
data. We find similar results for other events. Therefore, we 
conclude that there is no case of firm detection of additional 
planets. 


4.2. Constraints on Additional Planets 


Although no clear evidence of additional planets is found, 
it is still possible to constrain the existence of additional plan¬ 
ets around the lens. The constraint is expected to be strong 
because all analyzed events are highly magnified and thus ef¬ 
ficiency to planets would be high in wide ranges of planet 
parameters. 

For the purpose of constraining additional planets. 


2000 ; 


Albrow et al. 

2000; Kubas et al. 2008; Gould et al. 

2010; 

Cassan et al. 

2012h showing the confidence levels excluding 

the existence of an additional planet as a function of its sep- 


aration and mass ratio. To construct exclusion diagrams, we 
first estimate the detection efficiency of planets with respect 
to the planet parameters S 2 and qo according to the following 
procedure. 


(1) We first choose a triple-lens configuration with (s2 1 qi 1 '4>). 
For a given set of these parameters, we determine the re¬ 
maining parameters that best fit the observed light curve 
and then compute Ax 2 between the triple and binary fits. 


(2) We repeat the above process for many different orienta¬ 
tion angles %p of the second planet. 

(3) The detection efficiency for a given (S 2 ,q 2 ) is estimated as 
the fraction of the angles ip that produce significant devi¬ 
ations in the light curve. Then, the confidence level of ex¬ 
cluding an additional planet corresponds to the efficiency. 

As a criterion for planet detections, we set a threshold value 
A\ 2 h . iGould et~. ( '201 (f) suggested plausible values in the 
range A\ 2 h = 350-700 based on their experience in fitting 
lensing light curves for which planets were detected through 
the high-magnification channel. ICassan et all (120121) sug¬ 
gested a similar threshold. In our analyses, we adopt a median 
value A%i = 500. 

Figure [3] shows the exclusion diagrams of the individual 
events. For each diagram, the regions with different confi¬ 
dence levels of excluding additional planets are marked in dif¬ 
ferent colors. The star mark corresponds to the planet detected 
in the previous analysis. The upper x label represents the pro¬ 
jected planet-host separation in a physical unit (AU). The ticks 
marked by “J”, “S”, “U”, and “E” represent the planet masses 
corresponding to the Jupiter, Saturn, Uranus, and Earth of the 
Solar system, respectively. We note that the conversion from 
(si.qi) to the physical separation and mass of the planet is 
done based on the physical Einstein radius r E and total mass 
of the lens system M tot determined in the previous analyses, 
i.e., r±=sr E and M p . 2 = [qi/C + q\+qi)]Miot- 
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Fig. 3.— Exclusion diagrams showing the confidence levels excluding the existence of an additional planet with respect to the separation si and the mass ratio 
q 2 - For each panel, the upper x label represents the planet-host separation in physical units. The ticks marked by “J”, “S”, “U”, and “E” represent the planet 
masses corresponding to the Jupiter, Saturn, Uranus, and Earth of the Solar system. The star marks in the individual diagrams represent the planets detected from 
the previous analyses. 


In Table [2] we present the exclusion ranges of additional 
planets determinded with 90% confidence level for Jupiter, 
Saturn, and Uranus-mass planets. It is found that the con¬ 
straint on giant planets is considerably strong. On the other 
hand, the constraint on planets with masses less than Uranus 
is relatively weak. The weak constraint on low-mass planets 
is due to the small size of the planet-induced central caus¬ 
tic which decreases rapidly with the decrease of the planet 
mass dChung et al.l l2005h . In addition to planet masses, the 
sensitivity depends on various other factors. Some of these 
factors are related to observation, including the photometric 
precision, cadence, and completeness of the planetary signal 
coverage, etc. Other factors are related to the intrinsic prop¬ 
erties of planetary events, including the peak magnifications, 
severeness of finite-source effects, etc. 


For OGLE-2005-BLG-169, MOA-2009-BLG-319, and 
MOA-2010-BLG-477, the constraint is strong. These events 
have a common property of high peak magnifications with 
A max = 880, 208, and 599, respectively. 

By contrast, the constraint is relatively weak for 
OGLE-2005-BLG-071 (A max = 59) and MOA-2009-BLG-387 
(A max = 49) because the peak magnifications are low. For 
planets either with low masses on located away from the Ein¬ 
stein ring of the host, the central caustic is small and thus sig¬ 
nals of low-mass planets can be detected only for very high- 
magnification events where the source trajectories approach 
close to the caustic. 

For MOA-2007-BLG-400 and MOA-2008-BLG-310, the 
peak magnifications are high: A max = 629 and A max = 366, re¬ 
spectively. Nevertheless, the sensitivity is relatively low. We 
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find the main cause of the low sensitivity for these events is 
severe finite-source effects, which wash out planetary signals, 
especially for low-mass planets dBennett & Rhidll996l) . 

For MOA-201 l-BLG-293, on the other hand, the sensitiv¬ 
ity is not very high although the event was highly magnified 
(A max = 418) and experienced little finite-source effects. We 
find that the low sensitivity is mainly due to low photometric 
precision caused by the source faintness. 

5. CONCLUSION 

In order to investigate the potential existence of addi¬ 
tional planets, we reanalyzed high-magnification microlens- 
ing events for each of which a single planet had been previ¬ 
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ously detected. We found that introducing additional planets 
improves fits but the levels of improvement are not big enough 
to firmly identify additional planets. Although no clear evi¬ 
dence of additional planets was found, we could constrain the 
existence of additional planets in the parameter space. For this 
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